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ABSTRACT 

We present Spitzer 16 micron imaging of the Great Observatories Origins 
Deep Survey (GOODS) fields. We survey 150 square arcminutes in each of the 
two GOODS fields (North and South), to an average 3 sigma depth of 40 and 65 
/iJy respectively. We detect ~ 1300 sources in both fields combined. We validate 
the photometry using the 3-24 /xm spectral energy distribution of stars in the 
fields compared to Spitzer spectroscopic templates. Comparison with ISOCAM 
and AKARI observations in the same fields show reasonable agreement, though 
the uncertainties are large. We provide a catalog of photometry, with sources 
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cross correlated with available Spitzer, Chandra, and HST data. Galaxy number 
counts show good agreement with previous results from ISOCAM and AKARI, 
with improved uncertainties. We examine the 16 to 24 /im flux ratio and find 
that for most sources it lies within the expected locus for starbursts and infrared 
luminous galaxies. A color cut of S16/S24 > 1-4 selects mostly sources which 
lie at 1.1 < z < 1.6, where the 24 /im passband contains both the redshifted 
9.7 /im silicate absorption and the minimum between PAH emission peaks. We 
measure the integrated galaxy light of 16 /im sources, and find a lower limit on 
the galaxy contribution to the extragalactic background light at this wavelength 
to be 2.2 ± 0.2 nW m~ 2 sr" 1 . 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: 
high-redshift — infrared: galaxies 



Introduction 



UV light absorbed by dust is primarily reradiated in the far-IR (FIR), with a peak be- 
tween 60 and 100 /im. In addition, complex molecules, the polycyclic aromatic hydrocarbons 
(PAHs), radiate characteristic emission features in the mid-IR (M IR), the most promin ent 
of which are at wavelengths 6.2, 7.7, 8.6, 11.3 and 12.7 /im (see iPuget fc Legerl Il989l . for 
a review). Over the same wavelength range, there is continuum emission fr om very small 



dust grains, which can dominate PAH emission at wavelengths beyond 10 /im ( jLaurent et al. 
20001 ). The MIR flux, which results from the sum of these two emission mechanisms, cor 



relates strongly with the integrated IR luminosity from 8-100 /im, Ljr, which is a direct 
tracer of star formation (jKennicuttlll998l : IChary fc Elbazll200ll . and references therein). Ac- 
tive galactic nuclei (AGN) are also strong MIR sources, making the infrared an excellent 
tracer of obscured AGN which may not be accessible even to ultradeep X-ray observation s 



flLacv et al.ll2004i : ISaiina et al.ll2005l : IStern et al.ll2005l : iDaddi et al. 



Observations fro m the Infrared Spa ce Observatory (ISO; iKessler et al. 



adeep 
20071: 



Donlev et al.ll20Q8h . 



1996), the S pitzer 



Space Telescope (IWerner et al.ll2004l ). and the AKARI satellite ([Murakami et al.l 120071 ) have 
revolutionized the study of infrared luminous sources in the past decade. 

Of particular interest is the measurement of the integrated galaxy light and its value 
compared to the DI RBE measured extragalactic background light (EBL; for a review, see 
Lagache et al.ll2005l ). Lower limits on the EBL are inferred from the detection of galaxies 
in imaging surveys with, e.g. ISO and Spitzer (jElbaz et al.l l20Q2t iDole et al.ll2006l ). Upper 
limits on the background are obtained from 7-ray observations of blazars, because of the 
absorption of TeV photons by the cosmic infrared background through pair production (e.g. 
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Aharonian et al.l 120021 ) . More recently, iMatsuura et al.l ( 120101 ) have used long- wavelength 
(65-160 /xm) AKARI imaging to directly detect the cosmic infrared background, and confirm 
a value in excess of the lower limit measured by stacking of Spitzer sou rces. Considerable 



effort has gone into predicting the contribution of AGN to the EB L (IFadda et al.l 12002 



Barmby et al.l 120061 ; iTreister et al.l l2006t iBallantyne fc Papovichl 120071 . among others) . It is 
clear that galaxies hosting AGN do not dominate (well under 30%), but their precise con- 
tribution remains hard to estimate, in part due to the trouble identifying obscured sources. 
Measurement of the resolved portion of the background at a particular wavelength is per- 
formed by integrating the total flux from individually detected galaxies and through the 
technique of stacking. 

Deep, MIR observations of the Great Observatories Origins Deep Survey (GOODS; 
Giavalisco et al.l 12004 ; iDickinson et al.l 120031 ) fields have led to important advances in un- 
derstanding the global history of star formation and AGN evolution. For example, the 
peaks in the differ e ntial source count s at 15 and 24 microns at 0.4 and 0.2 mjy respectively 
felbaz et al1l2002[ Icharv et all I2004J : IPapovich et alll2004t iMarleau et aD l2004h. reflect the 
contribution of luminous infrared galaxies (LIRGs) at higher redshifts (e.g. iPearsonl 120051 ) as 
MIR features redshift into the bands. The GOODS fields also provide substantial informa- 
tion on the p resence of AGN, both obscured and uno bscured, using the ultradeep Chandra 
observations (IBrandt et al.ll200ll ; iGiacconi et al.l 120021 ). 



In this paper, we present Spitzer 16 /zm observations covering the GOODS fields. We 
provide the source catalog, including quality assessment flags, and provide source associa- 
tions between these observations and other available photometry from Spitzer, HST, and 
Chandra. We discuss preliminary analysis of the dataset including number counts, 16/24 
/jm colors, and the integrated galaxy light of 16 /im selected sources (including obscured 
AGN). Earlie r Spit zer imaging of a small region within GOODS-North was presented in 
Teplitz et all hoO$ ). 



Le Borgne et al.l (120091 ) used an earlier reduction of the present Spitzer survey as part 



of a study of the cosmic star-formation history using MIR number counts. They combine a 
non-parametric inversion of galaxy counts at 15-850 /im with constraints from measurements 
of the cosmic infrared background. They exclude a major contribution from "hyper-LIRGs" 
at high redshift, c onclu ding that these sources may in fact be AGN-dominated. In addition, 
Burgarella et al.l (120091 ) recently reported imaging of a portion of GOODS South at 15 jum 
using the Infrared Camera (IRC) onboard AKARI. Their analysis has some overlap with 
the Spitzer study (number counts, MIR colors), and the results are consistent, as we discuss 
below. While the two surveys are complementary, the Spitzer study had the advantage of 
significantly more telescope time; it covers more area and to greater depth, especially with 
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the addition GOODS-North. We report ~ 1300 objects over the two fields, compared to 
< 300 in the AKARI survey. 

We describe the survey, source extraction, quality assessment, and validation in Section 
2. We provide the catalog and discuss the survey properties in Section 3. In Section 4 we 
describe preliminary analysis of the data, before summarizing in Section 5. Throughout, 
we assume a A-dominated flat universe, with H = 71 km s" 1 Mpc" 1 , fi A = 0.73, and 
Q m = 0.27. 



2. Observations and Data Reduction 



In this section we describe the observations and data reduction. In addition, we have 
performed several validation checks and quality assessment measurements, which we describe 
in detail. 

The Spitzer peak-up imaging (PUI) capab ility offered a sens itive imaging capability at 
16 (Ltm using the Infrared Spectrometer (IRS; iHouck et al.l 120041 ) . The peak-up array was 
read out on the same detector as the Short-Low (5-14.5 /xm) spectroscopic channel, in a 
small field of view (56" x 80"). The Si:As detector was similar to that used in the 24 /im 
channel of the MIPS instrument (IRieke et al.ll2004l ). In 10 minutes of observation, the PUI 
achieved 5cr depths of ~ 45 /xJy 

Data were taken in two general observer Spitzer programs: GO-3661, observing GOODS- 
South in Cycle 1; and GO-20599, observing GOODS-North in Cycle 2. The southern survey 
was taken before the PUI mode was fully commissioned for Spitzer. In order to observe 
GOODS-South at 16 /im, the instrument was commanded to take a series of short spec- 
troscopic observations over a range of positions, and the images were acquired in parallel. 
Similar observations were taken over a small area within GOODS-North and are described 



in 



Teplitz et al.l (120051 ) . The northern survey used the standard PUI mode. 



The southern survey was designed to consist of ~ 130 square arcminutes of shallow data 
(2 minutes per pointing) and 10 square arcmin utes of deeper data (16 minutes per pointing) in 
the area of the Hubble Ultradeep Field (UDF; IBeckwith et al.ll2006l ). Observations were first 
taken in early 2005, but were compromised by persistent charge on the detector resulting from 
the preceding program, which targeted the rings of Jupiter. Much of the survey was repeated 
later in the year. Many exposures within the compromised observations were unaffected, so 
the shallow survey has an average of 4 minutes per pointing, with small areas of overlapping 
frames having greater depth. The UDF coverage varies from 16 to 32 minutes. 
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The northern survey covered 150 square arcminutes with ~ 10 minutes per pointing, 
observed i n 2006. We chose not to combine these data with the previous GOODS North 



imaging ( ITeplitz et al.l 120051 ) , given the small area of the latter and the different observing 
modes. 

IRS 16 /iin images were reduced by version S13.2 of the standard Spitzer Science Center 
(SSC) pipeline^. The pipeline supplied Basic Calibrated Data (BCD) frames with most 
instrumental effects corrected and flux calibration applied. Flux calibration was updated 
to match the latest version of the pipeline (S18.7). The pipeline removed a nominal low- 
background sky image, but some residual zodiacal light could have remained. We created 
median sky images from near-in-time subsets of the data, specifically from each contiguous 
block of observations or "astronomical observing request" (AOR). We then subtracted the 
median sky frames after scaling to the mode of the images. 

Individual PUI frames are quite small (56" by 80"), as the primary purpose of the 16 
/im camera was target acquisition. The plate scale is ~ 1.8 arcseconds per pixel. Geometric 
distortion is about 2%. 

We registered and combined images using the MOPEX software distributed by the 
SSC (IMakovoz fc Marleaull2005t iMakovoz et al.ll2006al Jbl). We employed drizzle interpolation 
( Fruchter &: Hook! 12002 ). MOPEX drizzle produces the same results as wdrizzle in IRAFfl 
The package uses the World Coordinate System (WCS) definition of both spatial offsets and 
geometric distortion. In imaging mode, the IRS did not perform an initial peak-up to refine 
the pointing, so the absolute WCS is good to only about 1 arcsecond, though the pointing is 
considerably better between BCDs within a single AOR. The Point Spread Function (PSF) 
at 16 /jm has a full width at half maximum (FWHM) of about 3.6 arcseconds. The final 
mosaics have a plate scale of 0.9"/pixel, and the pixfrac was set to 0.6. Figure [T] shows the 
mosaics of the GOODS fields. 

Photometry was performed by apply ing a custom poin t source extraction code which 
utilizes PSF-fitting and positional priors ( IChary et al.l 12004 ). We used 5a sources from the 
full GOODS IRAC channel 1 catalog (3.6 /xm; Dickinson et al. in preparation) as input. 
Sources within a 20" x 20" box around each input source were simultaneously fit to the 
measured PSF. The same technique was used in the preparation of the GOODS 24 /zm catalog 
(Chary et al., in preparation). The PSF for each field (North and South) was measured by 
registering about 15 bright point sources in each field. Source position was allowed to vary by 



1 http: / /ssc. spitzer. caltech.edu/irs/dh/ 

2 IRAF is distributed by NOAO, which is operated by AURA Inc., under contract to the NSF 
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up to an arcsecond to account for the uncertainty in the WCS. In addition to pointing issues, 
this flexibility has other advantages as well, including the possibility of actual displacements 
between the centroids of the starlight (IRAC) and the dust emission (MIPS or IRS), or cases 
where IRAC sources are blended (leading to a displacement of the centroid), but where one 
source may dominate the MIPS or IRS emission. We checked the residual map to ensure 
that no sources were missed by using the IRAC priors. 

The flux uncertainty was measured from residual pixels after subtracting the best-fit 
estimate of the source; the variance was taken as the PSF weighted sum of squares of the 
residuals. However, the noise in the drizzled images is correlated due to the sub-sampling. 
To c orrect for this effect, w e scaled the measured uncertainty upwards by a factor of ~ 1.7 
(see 



Casertano et al.ll2000l ). 



Photometric measurements are reported in Tables 1 and 3. We report the position of 
each object as that of the IRAC positional prior which was used as input to the source 
extraction. We cross correlated these positions with those of objects in the IRAC and MIPS 
24 /iin catalogs (Dickinson in prep, and Chary in prep., respectively) an d report their pho 



tometry as well. The IRAC photometry was measured with SExtractor (IBertin fc Arnouts 



19961 ). using aperture photometry except in the case of a small number of extended sources 



for which MAG _ AUTO is used. 

In general, we discarded sources with less than 5cr significance. This uniform cut allows 
us to consider the entire survey area, despite its non-uniform coverage. Final fluxes were 
consistent with aperture photometry, with appropriate aperture corrections. As described 
below, a small number of additional sources are reported in the catalog with less than 5cr 
significance (with a corresponding quality flag), because they were expected to meet that 
criterion based upon the integration time. 

Furthermore, to ensure a low incidence of spurious sources, we only used areas of the 
survey with good coverage. In the North and the UDF, we required at least ten individual 
data collection events (DCEs) per pixel. In the shallow Southern survey, this depth was 
not possible, but we did reject areas with less than 2 DCEs coverage. Depth of coverage is 
reported in the table. 

The IRAC observations of the GOODS fields are much deeper than the current survey. 
The 3.6 /iin catalog was used as input positional priors, but we can use the 8 jttm data as a 
quality check. Even accounting for the steep slope of IR luminous galaxies, it is unlikely that 
real sources will be detected at 16 //m and not at 8 jttm. We would reject as spurious any 
source which is not reported at 5a significance in the IRAC channel four catalog. However, 
all sources which meet our other selection criteria pass this test. Further quality assessment 
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is discussed below. 



2.1. Quality Assessment 

We performed several quality assessment procedures in order to flag objects which may 
be less reliable. Flags are reported in the table. 

First, we created a residual map by scaling the PSF to the measured fluxes and sub- 
tracting it from the image mosaic. We identified a small number of bright, extended objects 
(9 in the North and 7 in the South) for which the PSF-fit photometry is inappropriate. We 
measured these sources in a large aperture instead. A flag for extended sources is reported 
in the table. 

Second, we estimated the concentration of detected objects by measuring the ratio of 
aperture photometry in 6 and 2 pixel radii. Objects which are point sources should fall in 
the range of ~ 2 — 3 : 1. Objects outside the nominal range (ratio > 4 or < f .3) are flagged 
in the catalog. Several objects with anomalous values were inspected by eye and rejected 
from the catalog. 

Third, we flagged objects which may be affected by source confusion. We identified 
close companions in the catalog, within a radius of 5.4 arcseconds (6 pixels, a little larger 
than the FWHM of the PSF). We also flagged targets which have 8 /im sources (detected at 
> 3er) within 4 arcseconds, whether or not they are detected at 16 /im. 

Fourth, we compared the signal to noise ratio (SNR) estimate from our photometry 
code to the expectation from the exposure time map. The exposure time calculator at the 
Spitzer Science Center website predicts a la sensitivity of ~ 50 /iJy in a single 30 second 
exposure, and the observations are background limited so sensitivity is expected to scale 
with the square root of exposure time. The measured SNR generally agrees within about 
20% with expectations, though our estimates are slightly conservative. A small number of 
objects were selected as meeting the 5 sigma cut but exceeding the expected sensitivity; that 
is, they would have been excluded using the expected noise values. We retain these objects 
and flag them as long as the expected SNR is 4 or greater. Conversely, a somewhat larger 
number of objects were expected to meet the 5 sigma cut but did not in our measurements. 
We report such objects as low confidence sources as long as they have measured SNR greater 
than 4. 

Finally, we flagged sources with the spectral energy distribution (SED) expected for 
stars. Specifically, we identified stars as those objects for which the flux density is declining 
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across all six bands. We exclude one object for which / 2 4 > /s/2, which can be the case for 
some low redshift galaxies. Of these, 7 (2) in the N (S) are confirmed as stars in ground-based 
spectra. In the North, we also flag as stars two sources that are saturated in IRAC channel 
1 (which causes them to fail the SED check) but are confirmed in ground-based spectra. We 
check all star identifications in the ACS imaging of the GOODS fields and confirm that they 
are point sources. The imaging also shows that two of the objects identified as stars in the 
North are likely confused with nearby galaxies, and we flag these separately (see Section [3]). 



2.2. Photometric Calibration and Verification 

The calibration of the Spitzer Peak-Up Imaging mode is in the IRS Instrument Hand- 
book^] and briefly summarized here. The current calibration is based upon spectroscopic and 
imaging observation of four A stars for which detailed spectral models were available. There 
is a systematic uncertainty of ~ 5% in the calibration. 

Stellar fluxes were measured from the the IRS spectra of the calibrators by integrating 
under the filter transmission curve and color-correcting to the effective wavelength (15.8 
jum) assuming vF v = const. This assumption is designed to keep the color correction small 
for a wide range of spectral slopes. The calibration is normalized to infinite aperture using 
Tiny Tim V2.CQ In practice, aperture photometry requires an aperture correction. Profile- 
weighted fitting, such as that utilized in the present study, is normalized to a finite radius 
and then aperture corrected as well. 

The PUI calibration is tied to the calibration of the spectrometer, which has been 
validated against a large number of stars. To evaluate the calibration in the particular case 
of the GOODS data, we can compare the photometry of stars within the two fields. Twenty- 
five objects have IRAC colors clearly indicative of stars. However, four of these objects have 
quality flags greater than 1, so we exclude them. We also exclude another 6 objects that 
have SEDs indicating possible IR excesses. So, we have 15 stars with which to validate the 
PUI calibration. These objects are faint, however, and we do not know their stellar types. 
Nonetheless, the mid-infrared spectral energy distribution of stars does not vary much across 
a wide range of stellar types. The Spitzer Atlas of Stellar Spectra (SASS, Ardila et al. in 
preparation) observed stars across many stellar types with the IRS. In Figure [2] we plot the 
photometry for the 15 stars compared to SASS IRS spectra of stars with types A through 



3 http://ssc. spitzer. caltech.edu/irs/irsinstrumenthandbook/ 

4 http://ssc. spitzer. caltech.edu/archanaly/contributed/stinytim/index. html 
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M, excluding super giants, and the iKuruczl (119791 ) model for an AOV star. We apply a 
small (< 4%) color-correction to the photometry to account for the difference between the 
assumed calibration reference spectrum [yf v = const for IRAC0 and PUI; 10 4 K black body 
for MIPS0) and a typical (5000 K black body) stellar spectrum. The spread in photometry 
is similar to the expected variation within stellar types, and does not indicate systematic 
offsets in the photometric calibration. 



2.3. Comparison with ISOCAM and AKARI 



In addition to validation using stars, we compare Spitzer photometry to other observa- 
tions of the field at similar wavelengths. 



ISOCAM (jCesarsky et al.l If 9961 ) on-board the Infrared Space O bservatory observed a 
portion of GOODS-North centered on the Hubb le Deep Field North ([Williams et al.lll996l ). 
Careful reduction of those data was presented by lAussel et al.l (119991 ) . The ISOCAM catalog 
reports 40 objects detected at 15 /mi in the main catalog. Of these, 36 are associated with 
sources in the Spitzer 16 /im catalog. The four unmatched sources are all fainter than 100 
/iJy, and two of them do not appear in the (deeper) 24 am catalog, indicating that they 
are probably spurious. Figure [3] shows the comparison of Spitzer and ISOCAM fluxes for 
matched sources. 

Both ISOCAM and Spitzer PUI are calibrated assuming vf u = const, but the effective 
wavelengths are substantially different, being 14.3 and 15.8 /im, respectively. So, a color- 
correction to the ISOCAM photometry equal to the ratio of the effective wavelengths - e.g., 
a factor of 1.1 - is expected. We fit the offset between Spitzer and ISOCAM photometry 
assuming a constant slope of unity and find Spitzer fluxes to be 1.13 times brighter. Overall, 
we see general agreement within the large error bars. Spitzer had better spatial resolution 
at 16 /im than ISOCAM, and some of the ISOCAM sources are blends of multiple objects, 
and this may explain the difference in flux for some objects. The four brightest objects do 
not appear to suffer from blending, however. 



As noted in lTeplitz et al.l ( 120051 ). the difference in filter bandpasses can cause substantial 
differences in the reported flux densities at some redshifts (up to a factor of 1.5 or even 2) as 
prominent MIR features move in and out of the filters. Figure H] shows the ratio of Spitzer 



5 IRAC Instrument Handbook; http://ssc.spitzer.caltech.edu/irac/iracinstrumenthandbook/ 
6 MIPS Instrument Handbook; http://ssc.spitzer.caltech.edu/mips/mipsinstrumenthandbook/ 



7 See the "ISOCAM Photometry Report" , 1998; http://www.iso.vilspa.esa.es/users/expldib/CAM/photom_rep_fn.ps.gz 
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to ISOCAM fluxes for objects with known redshifts (see Section [3T2T) . 

The four brightest objects in common between Spitzer and ISOCAM show a larger offset 
in photometry than is typical for the sample as a whole. The differences may result, in part, 
from the difference in filter bandpasses as shown in Figure HI All four objects are detected 
by Chandra (see Section [373]) . two of them in the hard band. The sources do not appear to 
be variable, howeve r, as all four were also detected by Spitzer in independent measurements 



( iTeplitz et al.l 120051 ) . and their flux densities agree within 10%. 



Burgarella et al.l ( 120091 ) used the IRC on-board the AKARI satellite to observe a ~ 
10' x 10' region partially overlapping with GOODS-South. Their catalog contains 67 objects 
detected at > 5a in areas which we observed with sufficient coverage (at least 2 DCEs). 
Of these, 60 have associated sources within 3" (two IRC pixels) in the Spitzer catalog. 
Figure shows the comparison of Spitzer and AKARI flux densities for matched sources. 
AKARI IRC is also calibrated assuming vj v = const, but with an effective wavelength of 
15.0 jum, implying a color correction between the filters of 1.06. However, the best fit to the 
photometry shows a difference in photometry of a factor of 1.3. The filter bandpasses are 
slightly different, with the IRC filter being considerably wider, which may account for some 
of the difference. 



3. Results 

Extracted sources are given in Tables 1,2 and 3,4 for GOODS-N and GOODS-S, respec- 
tively. The tables include ancillary data from other telescopes, together with other Spitzer 
photometry in the IRAC and MIPS 24 jum channels. Quality flags are given in the tables; 
caution is recommended in the interpretation of results based upon sources that do not have 
quality 1. 

Specifically, the tables provide the available data on each source: position, photometry 
and uncertainty from Spitzer, HST, and Chandra, spectroscopic redshifts where available, 
and quality flags. The process for identifying HST, Chandra and redshift associations is 
given later in this section. Details of the columns in the table are given below: 

• Column (1) gives the source ID number within the catalog. 

• Columns (2) and (3) give the right ascension and declination (J2000) of the source. 
Positions are reported as the IRAC prior position used as input to the source extraction. 

• Column (4) gives the coverage (in number of exposures) of the central pixel of each 
source. Exposure times were 30 (N) and 60 (S) seconds. 
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• Columns (5) through (16) give the IRAC channels 1-4, IRS 16 /an, and MIPS 24 /im 
flux densities and uncertainties for the source, in units of /xJy. Uncertainties are given 
after each photometric point (e.g. Column 5 is photometry, Column 6 is uncertainty). 
The IRAC uncertainties include a 5% systematic uncertainty in the calibration (added 
in quadrature to the measured uncertainty). This systematic term dominates the 
uncertainty for most objects, especially in channels 1 and 2. 

• Column (17) gives the spectroscopic redshift for the source, where available. Spectro- 
scopic counterparts were usually chosen to be the closest optical/NIR object to the 
16um position. In a few cases, the 16um emission is likely due to a blend of sources, 
and the slightly more distant one is the more likely dominant counterpart. A value of 
-99 indicates that no redshift is available. The very few stars with spectra are listed 
with z=0. 

• Column (18) gives the reference code for the spectroscopic redshift (see Section l3T2~j) . 
Redshifts were collated from the literature, as well as from observations by the GOODS 
team (Stern et al. in prep.). 

• Column (19) indicates that the source has non-zero flux in at least one Chandra band 
(a value of 1 is detected). 

• Column (20) and (21) give the soft and hard X-ray fluxes, respectively, for the associ- 
ated Chandra detection. 

• Column (22) is a star flag, with a value of 1 indicating the source is likely a star. A 
value of 2 indicates that the position of the IRAC prior is confused between a star and 
a neighboring galaxy (there are two such cases in the Northern field). 

• Column (23) gives the 16 /im concentration index, defined to be the ratio of flux within 
apertures of radii 6 and 2 pixels. 

• Column (24) gives the number of 8 /im sources within 4 arcseconds. 

• Column (25) gives the number of 16 /tm sources within 5.4 arcseconds. 

• Column (26) gives the bit-wise quality flag for the source: Bit - the object was 
determined to be real and included in the catalog; if only this bit is set, then there are 
no notes and the object has the best quality; Bit 1 - possible confusion because more 
than one 16 /j,m source lies within 5.4 arcseconds; Bit 2 - bad concentration index; Bit 
3 - low SNR but the source was included because the expected SNR was at least five. 
In GOODS-S, we also include a flag for Bit 4, in the case of sources with coverage of 
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two exposures instead of three. Because flags are assigned bit-wise, sources may have 
multiple flags set. So, for example, an object with a flag value of 19 would indicate that 
Bits 0,1, and 4 had been set because: (Bit 0) the object was included in the catalog; 
(Bit 1) the measurement of the object may suffer from confusion due to a close 16 fim 
neighbor; and (Bit 4) the object had a coverage of only two exposures. 

• Column (27) flags extended sources which were measured with aperture photometry 
instead of the PSF-fit. A value of 1 indicates that the source was extended. 

• Column (28)-(31) give the HST magnitudes for the B,V,I,z bands respectively (see 
Section 13. lj ). 

• Column (32)-(35) give the HST uncertainty (in magnitudes) for the B,V,I,z bands 
respectively. 

• Column (36) gives the number of I-band sources within 1 arcsecond radius. 
We report the detection of 840 (North) and 476 (South) sources. 

The depth achieved in the survey is consistent with expectations. The 5a limit in the 
shallow Southern area varies from ~ 65 to 85 /xJy. The UDF limit is ~ 30 /iJy, though the 
small area means that there are few objects even in the faintest bin. The Northern survey 
reaches ~ 40 /iJy. Figure |6] shows the area covered in bins of predicted sensitivity. Table 
[5] summarizes the number of objects in each field, and in some sub-categories (see below). 
Table |6] lists the predicted and achieved depths for the North, South, and UDF regions. 

In the remainder of this section, we examine the association between 16 /im sources and 
other measurements within the GOODS fields. 



3.1. Optical Photometry 



We use the optica l photometry from the publicly released GOODS version 2.0 catalogs 
( IGiavalisco et al.ll2004[ )l 8 l. We use the reported MAG_AUTO fluxes. We cross correlate the 
positions of 16 /im detections with those of optical sources within 1 arcsecond. The GOODS 
catalogs used the z-band images for source detection and then measured photometry in the 
other three optical bands. 



http: / /archive. stsci.edu/prepds /goods / 
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Of the 840 sources in the north, 809 have associated optical sources within 1 arcsecond. 
In the South, 465 of 476 sources in the 16 /im catalog have optical associations. Of the 31 
(N) and 11 (S) sources without an associated optical detection, 11 (N) and 4 (S) do not 
have coverage in the ACS mosaics. For objects with coverage, 6 (N) and 4 (S) are near the 
diffraction spikes from bright stars, and 6 (N) and 1 (S) are in the outskirts of extended 
galaxies. The remaining 8 (N) and 2 (S) sources may be very red, or they may have a larger 
uncertainty on the centroid - either due to Spitzer pointing uncertainty or to blending with 
another source in the IRAC positional priors. For matched sources, the possibility of source 
confusion is significant when matching MIR and optical sources. We find 57 (N) and 37 (S) 
sources have multiple possible optical matches within 1 arcsecond, and we report the closest 
one to the IRAC position. The number of close matches is flagged in the table. 



3.2. Redshift Distribution 

Spectroscopic redshifts have been measured for sources in the GOODS fields by numer- 
ous surveys. The specific spectroscopic surveys used are given in the tables. 

We correlate the positions of detected Spitzer sources with redshift identifications within 
about 1 arcsecond radii. In cases for which multiple optical sources are within 1 arcsecond 
of the 16 /im target, we usually take the closest positional match. In a few cases the 16 /im 
emission is likely due to a blend of sources, and the slightly more distant one is the more likely 
dominant counterpart. Another 21 (N) and 8 (S) objects have optical sources in the redshift 
catalogs between 1 and 2 arcseconds away, allowing the possibility for misidentification in 
those well. Redshifts are reported in the table. In the North 701 redshifts are 

available out of 826 objects not identified as, or confused with, stars. Most (654) of the 
redshifts are for objects with %ab < 24, and at that magnitude, redshifts are available for 
94% of 16 jum sources. In the South, 381 redshifts are available out of 466 galaxies detected. 
Only 36 of the spectroscopic redshifts in the South are for objects fainter than %ab = 24, 
and 87% of 16 /im targets down to that magnitude have redshifts. For comparison, there 
are 828 redshifts associated with 5a 24 /im detections in the North and 814 in the South. 

Figure [7] shows the redshift distribution of detected 16 /im sources. It appears that the 
distribution of associated redshifts is similar for 24 and 16 /im sources at redshifts less than 
1.5. At higher redshifts, there are significantly more 24 /im detections, as a result of both the 
greater sensitivity of the MIPS observations and the intrinsic brightening of starburst sources 
when the strong PAH features shift into the 24 /im band. The similarity of the distributions 
at lower redshifts suggests that the distribution may be a stronger function of the optical 
limits on obtaining redshifts rather than the MIR observations themselves. Figure [8] shows 
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the 16 /im flux density versus redshift. 

The median spectroscopic redshift, excluding stars, for 16 fim sources is 0.85 (N) and 
0.82 (S), and the mean is 0.86 in each field. About 30% of redshifts in the North are at z > 1, 
and about 35% in the South. About 2% of the sources in each field (13 North, 6 South) have 
z spec > 2.0. Most of these (~half in the North and 5/6 in the South) are identified as AGN, 
as described below. The highest spectroscopic redshifts for 16 /im sources are 3.48 (N) and 
3.47 (S). 



3.3. X-ray Sources 



The Chandra 2 Msec surveys of the GOODS North and South fields, respectively, are 
the deepest X-ray observations tak en to date. We c o mpar e the 16 fim survey to the Chandra 
catalogs of klexander et all J200J and ILuo et all hoO$ ). There are 308 (N) and 293 (S) 
X-ray sources within the GOODS area surveyed at 16 /xm. Of these, 117 (N) and 92 (S) 
are associated with 16 /zm sources, comprising 14% and 20% of the 16 /im samples. Sources 
that are detected in the Chandra catalogs are flagged in the table. 



4. Discussion 

4.1. Active Galactic Nuclei 

The detection of sources at 16 fim selects both strongly star-forming galaxies and those 
which host AGN. Most sources in the present survey lie at redshifts below 2, which ensures 
that the PUI passband samples the wavelength range that covers emission by dust rather 
than direct stellar light. At z > 1, the PUI band is increasingly dominated by hot dust 
characteristic of AGN, and at z > 2 the sensitivity of the survey is mostly limited to such 
sources. 

Hard X-rays are usually indicative of an AGN. The ultradeep Chandra data select many 
of the AGN within GOODS. There are 96 (N) and 58 (S) PUI sources associated with 
detections in the hard band. Given the depth of the CDF s, a few of th e se ob jects may 



have X-rays from purely star-forming galaxies. For example, iTreister et all (120061 ) excluded 
sources with Lx < 10 42 ergs s _1 when selecting AGN in GOODS. The fraction of such 
sources is small amongst 16 /im targets, and the exact nature of those sources is uncertain, 
so we retain them as possible AGN in this analysis. 



Some galaxies hosting AGN, even those whose bolometric luminosity is dominated by 
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them, are undetected even by Chandra. These sources may be selectable using their MIR 
color, when AGN heat ed dust dominates over stellar light. In many cases, these will be have 
power-law SEDs (e.g. iDonley et al.l 120071 ). We perform a x 2 fit to the IRAC+16 /im colors 
to identify power-law sources. We fit all SEDs with a power-law, u~ a , and select those with 
a > 0.25. We do not set a limit on the goodness-of-fit, but we exclude sources where the 
channel 3 or 4 flux does not exceed the channel 1 flux. Using only IRAC photometry with 16 
(ttm does not significantly change the results. We find about 10% more AGN when combining 
X-ray with power-law selection than when using x-ray selection alone. Conversely, we find 
that about 45% of hard X-ray sources have power-law SEDs. 



A similar fraction was measured by iBarmby et al.l (120061 ) . who found that ~ 40% of 
X-ray sources in the extended Groth Strip have red power-law SEDs in the IRAC bands. 

Combining the two selections, we find 105 (N) and 65 (S) AGN. These numbers are 
likely to be a slight underestimate, as they may not include Type 2 AGN whose IRAC colors 
are not true power-laws.. Nonetheless, we find that about 15% of 16 /zm sources are galaxies 
hosting AGN. Figure shows the fraction of sources with AGN in bins of 16 /im flux density. 
The fraction is higher for the brighter sources. 



Amongst rare, bright sources, obscured AGN may be more common. IStern et al.l (120051 ) 
identify QSOs and Seyfert 1 galaxies as objects with red IRAC colors, but their selection 
is conta minated by star for ming galaxies at the faint fluxes, such as those in the GOODS 



Donley et al.l (120071 ) find that when using the IRAC color-color selection only 55% 



survey. 

of MIR power-law AGN are detected in the X-ray. 



Treister et al.l ( 120061 ) find a decreasing AGN fraction with decreasing 24 /im flux for 
sources in the GOODS fields, down to ~ 8% at < 100 /iJy using purely hard X-ray selection. 
The fraction of X-ra y sources at 16 /im i s similar, as expected given the high rate of detection 
of Chandra sources. IDonley et al.l ( 120081 ) also find a decreasing AGN fraction with decreasing 
24 /im flux, down to about 10% X-ray AGN. They then expand the estimate to MIR-selected 
AGN and estimate about 15-25% at 100-300 /iJy, similar to our estimate given the large 
error bars. 



Table [5] summarizes the numbers of AGN per field. 



4.2. 16 to 24 urn Color 



One of the most prominent MIR spectral features is the broad (full width ~ 2 /mi) 
silicate absorption trough at 9.7 /mi. Attenuation at this wavelength can approach an 
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order of magnitude in typical ultraluminous IR galaxies (e.g. lArmus et al.l 120041 ) . As a 
resu lt, this feature will signifi cantly depress the photometry measured in a broad-band fil- 



ter. iTakagi &: Pearson! (120051 ) suggested that the 16 to 24 fim ratio can be used to identify 
"silicate-break" galaxies, at redshifts 1.1 < z < 1.6 where the silicate absorption is solidly 
within the MIPS 24 /xm bandpass. 

Figure [10] shows the S^/S^ ratio of catalog sources for which spectroscopic redshifts 
are available, and which are solidly detected in the 24 /iin band. We apply a small color 
correction (4%) to the MIPS 24 /iin flux densities in order to account for the difference in the 
way the instruments are calibrated. The PUI is calibrated assuming a reference spectrum of 
vF v = const, but MIPS is calibrated using a 10 4 K black body (see the MIPS Instrument 
Handbook). 

Much of the scatter in the flux ratio is expected due to variation in the source SEDs. 
We demonst r ate th is by calculating the expected ratio for redshifted local template sources. 
Smith et al.l ( 120071 ) obtained Spitzer IRS spectra for a range of local starbursts, and found 
a large variation in both PAH emission strength and silicate absorption depth. They find 
that the relative strength of features can vary by a factor of two in star forming galaxies 
(Lir < 3 x 10 11 L Q ) depending upon source propert ies, though they no te the variation 



should be somewhat less in more luminous sources. iBrandl et al.l (120061 ) also measured 
local starbursts, notably NGC 7714, and calculated an average starburst spectrum, which 
is redder than the Smith et al. templates, and so results in a lower expected ratio. In the 
plot, we indicate the expected ratio for these template s as a shaded region. Similar results 
have already been reported by iBurgarella et al.l (120091 ) using the AKARI measurements in 
GOODS-S. They find that local starburst templates mostly reproduce the flux ratio well, 
but that there are some discrepancies at higher luminosities. The AKARI survey was limited 
to z < 1.2. 



Armus et al.l (120071 ) calculated the expected 16 to 24 am ratio for local ultraluminous 



IR galaxies (ULIRGs), and found that it is significantly lower than for starbursts at most 
redshifts because ULIRGs have a red continuum. However, in the redshift range for which 
the silicate absorption falls in the 24 /iin passband, the ratio can be extreme (> 3 for some 
ULIRGs). In the Figure we plot the ratio using Arp 220 as a representative ULIRG template. 
Finally, we note that AGN lacking strong silicate or PAH features will have red colors at all 
redshifts. Figure [10] also shows the ratio for an AGN with a very red (a = —2) power law 
SED (/„ oc v a ). 

For clarity, error bars are omitted from points in the figure within the expected region. 
Typical uncertainties are 0.1-0.2 in the flux ratio. Instead, we only plot individual uncer- 
tainties for objects with 16/24 colors more than la bluer than the expected regions. There 



-17- 



is significant scatter in the plot, but a trend at z > 1 indicating the presence of silicate 
absorption is present. 

A number of sources at z < 1.1 are anomalously blue, though most of them are likely the 
result of measurement error (large uncertainties) or data quality issues (most have quality 
flag warnings, as shown in the figure). The rest may have SEDs that are bluer than local 
starbursts because the 24 /im band is tracing longward of 12 /im emission and their warm 
dust contribution is low. 

Objects at 1.1 < z < 1.6 with a ratio bluer than .SW ^ = 1-2 lie solidly in the 



iZ!6 

region expected for silicate absorption. iKasliwal et all (120051 ) used the same selection when 



identifying silicate absorbing sources in early Spitzer observations of the NOAO Deep Wide 
Survey field. However, many objects at lower redshifts fall within the same cut. Selecting 
objects with a ratio > 1.4 eliminates many, though not all, of the low redshift interlopers. 
If we consider only sources with a 16 /im quality flag of 1 and an available spectroscopic 
redshift, then 10 of 44 selected sources are outside the redshift range for silicate absorption. 
For comparison, ~80% of 16 /im sources with spectroscopic redshifts are at z < 1.1. If 
we consider all sources which are securely detected at 24 /im and have available redshifts 
(regardless of quality flag), then 33 of 91 sources are possible interlopers; however, about 
half of these have quality flag with either Bit 1 or Bit 2 set, indicating possible confusion 
or a bad concentration index. We note that not all objects with blue ratios are true silicate 
absorbers. Many, perhaps most, are star-forming galaxies wher e the blue color res ults from 



the dip between the 7.7/8.6 and 11.2/12.7 /im PAH complexes (ISmith et all 120071 ) 



A few objects lie at redshifts higher than the range where silicate absorption falls in 
the 24 /im band but still have colors bluer than expected for either AGN or the starburst 
template. These objects could have buried AGN causing large hot dust emission at 5-6 /im, 
or they could have anomalously strong 6.2 PAH features. 

Figure [11] shows the flux ratio for objects selected as AGN by either power-law SEDs 
or hard X-ray detection. The S\§/S2A > 1.4 selection does not appear to select objects that 
are dominated by known AGN. In the Figure, we also see a number of AGN sources with 
very low S^/ S24 ratios. At z ~ 0.6, some of these may be explained by silicate absorption 
in the 16 /im band. Of course, the ratio will be low for objects with red SED as well, so the 
silicate selection is not as clean. 

Figure [T21 shows histograms of the flux ratio for the 16 /im sample, and the AGN subset. 
Applying the ratio cut at 1.4 (and excluding sources with other spectroscopic redshifts), we 
find 107 blue sources, of which possible 68 have a quality flag of 1, which could indicate that 
they are z ~ 1.3 sources. Table [5] compares these numbers with the AGN selected for each 
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field. 



4.3. Number Counts 

We calculate galaxy number counts for each field. First, we determine the area, Ai, over 
which each source could have been detected assuming our estimate of the sensitivity based 
upon the depth of coverage. Then, we sum the reciprocal of the area for each source within 
logarithmically spaced bins of flux and divide by the binwidth, 5f. 

nc = 77T\ I^T" W 

{OJ) bin f-** Ai 

This calculation does not include the effects of source confusion which could lead to 
slightly under-counting the number of galaxies. We estimate the uncertainty on the counts 
for each bin as the Poisson error in the measurement of the number of sources in the bin. 

Comple t eness corrections were calculated using a Monte-Carlo simulation, following 



Chary et all ( 120041 ) . Separate simulations were performed for the three depth tiers of the 
survey - GOODS-North, GOODS-South, and the UDF. Artificial sources were added to the 
original data images and recovered. These sources were placed at random positions. The 
flux distribution of simulated sources was flat in log(f u ). The fluxes varied from 20 to 1000 
fjjjy. To avoid confusion of simulated sources, only a small number were added at a time and 
the simulation was repeated many times. In the North and South, 50 objects were added to 
the images at a time, with the simulation repeated to build up ten thousand simulated input 
sources. In the UDF, only 15 sources were added at a time. Sources were recovered using 
the same positional prior code used for the catalog generation, discarding sources below 5a. 
The completeness of recovered sources is seen to be worst at the faint fluxes, as expected. 
There is considerable incompleteness at relatively bright fluxes, as well, due to confusion of 
simulated and real sources. This effect appears to be on the order of 5-10%. 

A matrix Pij for the output flux distributi on of the artificial sources was generated, where 



is the input flux and j is the recovered flux ( Smail et al.lll995l ). The classical completeness 



i 

for sources in the i th bin is the ratio of number of recovered sources in that bin to the sum 
over all j for that bin. The observed catalog of sources in the real image was then distributed 
among the flux bins. The P{j matrix was renormalized such that the sum over i for each j 
was equal to the number of detected sources in that flux bin. The completeness corrected 
counts in each flux bin i is then the sum over j of the renormalized Py matrix. 

In the North, we estimate completeness to be ~ 50% at 40 /iJy (the approximate limit 
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of reported sources), and 80% at 60 /xJy. In the South, we estimate 80% completeness at 
80 /iJy in the shallow survey, with a steep fall off at fainter fluxes due to poor coverage. In 
the calculation of the number counts, we do not use the shallow Southern survey below 65 
/xJy. In the UDF, we estimate 50% completeness at 30 /iJy (the limit of reported sources) 
and 80% at 50 /iJy. 

Figure [13] shows the measured number counts, including the completeness correction 
described below. The figure al so compares 15 u r n num ber counts from ISOCAM and AKARI, 



but without color-correction. iLe Borgne et al.l (120091 ) showed that these number counts are 



dominated by low redshift objects (z < 0.5) at fluxes above 200 /iJy, and by moderate 
redshift ones (0.5 < z < 1.5) at fainter levels; they infer a small contribution (< 20%) from 
higher redshift sources in bins below 100 /iJy. 

In the figure, we see a significant difference between the counts in the Northern and 
Southern fields. The counts peak around 0.4 mJy in the North, but around 0.2 mJy in the 
South. It is likely that this effect is the result of cosmic variance. However, the distribution 
of redshifts (see Section 13. 2p sh ows no clear ev i dence for an over-density in the Northern 



field to explain the difference. iLeFloc'h et al.l (120091 ) noted a similar effect in comparing 



GOODS-North counts at 24 /zm to those in the (wider area) COSMOS field. 



Elbaz et al.l (119991 ) measured a faint-end slope of a = —1.6 from the ISOCAM 15 /im 
counts. They defined the faint-end to be bins with S15 < 0.4 mJy. In the same range, we 
measure a slightly steeper slope of -1.9 when fitting the counts from both the North and 
South. However, the fit is dominated by the high significance bins between 0.1 and 0.4 mJy, 
where the two fields have different peaks. If we consider only S\q < 0.2 mJy, then a = —1.7 
is consistent with the data. 

In Figure [HJ we show the contribution of AGN to the number counts. The result is 
largely similar to Figure [91 with X-ray selected AGN dominating the AGN contribution at 
bright fluxes, and the total AGN contribution decreasing at faint fluxes. We also show in 
the figure the contribution of silicate absorption candidate sources, which could be highly 
obscured AGN. These objects occur mostly at faint flux levels, given their high redshift. 



4.4. Integrated Galaxy Light at 16 microns 



Elbaz et al. 


2002; 


Altieri et al. 


1999 


)• 



contribution of ISO-detected galaxies to the 15 /im EBL to be 2.7 ± 0.62 nW m sr at 
S*i5 > 30 /iJy by integrating the flux measured from faint sources (the integrated galaxy 
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light, IGL) including the observations of lensing clusters. iHopwood et al.l (120101 ) improved 
the measurement of the lensing cluster and determined an IGL value of 1.9 ± 0.5 nW m -2 
sr -1 at Sis > 10 /iJy. 

We can improve the measurement of the ~ 15 /im IGL using the Spitzer sources. We 
have significantly more detections fainter than 50 /iJy than any of the ISOCAM surveys, and 
we cover two independent fields that were covered by the ultradeep AKARI measurement. 
We first take the average of the number c ounts for both GOODS fields, and combi ne them 
with ISO CAM counts at flu xes > 1 mJy felbaz et al.lll999l : iGruppioni et al.l l2002h . Next, 
following lElbaz et al.l ( 120021 ). we integrate dIGL/dS (defined by their Equation 6) over the 
range 30 /iJy to 1 mJy, after fitting the counts with a 3 rd degree polynomial (the dotted line 
in Figure IT3"j) . A conservative estimate of the uncertainty on IGLi 5 was obtained by fitting 
and integrating the Poissonian upper and lower la uncertainty on the counts. 

We find IGL W (S w > 30/iJy) = 2.2 ± 0.2 nW m~ 2 sr^ 1 . Note that if we color-correct 
the Spitzer sources to the effective wavelength (14.3 /an) of the ISOCAM LW3 filter as in 
Figure [31 we would obtain I GL^ (Sy\ > 25 / iJy) = 1.8 ± 0.2 nW m~ 2 sr -1 , just outside 



the la uncertainties from the iMetcalfe et al.l ( 120031 ) result. Figure [J5] shows the IGL as a 
fun ction of sensitiyity li mit. The upper limit on the EBL, 4.7 nW m -2 sr -1 is taken from 
the lRenault et al.l (120011 ) high energy 7-ray measurement of Mrk 501. Comparing the upper 
limit to the EBL, following Metcalfe et al., we find the 16 /im counts down to 30 /iJy appear 
to be resolving at least ~ 50% of the monochromatic 16 /im EBL. Most likely the resolved 
fraction of the 15 /iEBL is even larger since the IGL at the faintest flux limits probed by 
the Spitzer surveys (Figure 15) appears to be asymptoting towards a value of ~3 nW m~ 2 
if extrapolated to zero flux. 

In principle, the 16 and 24 /im IGL should be similar as they are sampling mostly similar 
populations. Though, the 24 /im band is more sensitive to higher redshifts and, as we have 
shown, the 16 /mi band picks up a slightly higher fraction of AGN . Our value is lower than 



the quoted estimate of the 24 /im IGL obtained by lPapovich et al.l (120041 ) of 2. nW m 2 



sr . However, that result was obtained by measuring the 24 /im number counts down to 60 
/iJy and extrapolating to fainter fluxes using a fit to the faint-end slope. If we extrapolate 
the 16 /im source counts to fluxes 10 times fainter (down to 3 /iJy), assuming faint-end slope 
1.9, we obtain IGLiq ~ 2.9 nW m~ 2 sr -1 , in good agreement with the MIPS value. 

Fadda et al.l (120021 ) find that AGN contribute ~ 17% of the 15 /im background (see 



a 



also lLa Franca et al.l 120071 ) . Similarly, iMatute et al.l (120061 ) place a lower limit of the AGN 
contribution to the 15 /im IGL of 4-10%, using only optically-selected AGN. We measure 
the IGL for the AGN among the 16 /im GOODS sources, and find a contribution to the 16 
/im IGL of ~ 15% from X-ray and power-law selected AGN. 



-21 - 



Spitzer measurements at 8 an d 24 /jm have typically d eterm ined a marginally smaller 
contribution to the IGL from AGN. iBallantyne fc Papovichl (120071 ) and lBarmby et al.1 (120061 ) 
find contribution of ~10% to the 24 /im background base d on X-ray sele c tion o f AGN. Most 
of these surveys also rely on X-ray selection of AGN, and iBarmby et al. suggest that 

the AGN con tribution be c orrec ted upwards by a factor of 1.5 to account for Compton- 
thick sources. Brand et al.1 (120061) s ugges t that only 3-7% of the 24 /im background results 
from AG N. and iFranceschini et al.l ( 120051 ) place the contribution at 10-15% based on optical 
selection. ISilva et al.l (120041 ) note that the fractional AGN contribution is, itself, an upper 
limit to the contribution to the IGL from AGN radiation because significant IR light from 
these sources arises from the host galaxy. 



5. Summary 

We have presented source catalogs for Spitzer 16 /im observations of the GOODS fields 
using the Peak-up Imaging capability. We surveyed the ACS area (150 sq. arcminutes) to 
depths of 40 and 65 /xJy (50% completeness) in the Northern and Southern fields, respectively. 
In the ~ 10 sq. arcmin of the UDF, we reach 30 /iJy. We detect 840 (N) and 476 (S) objects. 
These observations are the widest contiguous-area PUI taken during the Spitzer mission, and 
among the deepest. The catalog will be available through the NASA/IPAC Infrared Science 
Archive (IRSA@). 

We validate the photometry by demonstrating that the Spitzer SEDs of stars in the 
fields are consistent with standard templates from Spitzer spectroscopy. We compare PUI 
photometry with observations within GOODS North by ISOCAM and find reasonable agree- 
ment after color correction, though the uncertainties are large. Comparison with AKARI 
shows marginally significant disagreement (~ 30%), with Spitzer fluxes being higher. We 
match PUI sources within 1 arcsecond with Spitzer, Chandra and HST detections. We also 
report spectroscopic redshifts from the literature where available. The distribution of red- 
shifts is similar to that for 24 /tm sources despite the shallower MIR depth, suggesting that 
redshifts so far available are mostly limited by optical faintness. We flag sources in the cat- 
alog that may have quality issues: possible source confusion, bad concentration index, lower 
than expected signal-to-noise ratio, low coverage. 

Preliminary analysis of the survey find the following: 
• Matching 16 fim sources with Chandra hard-band detections finds that ~ 11% of PUI 



http: / /irsa. ipac.caltech.edu 
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sources have X-ray counterparts. Combining these sources with power-law MIR selec- 
tion for obscured AGN finds that about 15% of objects in the catalog are potentially 
AGN. The fraction of AGN increases with increasing 16 /im flux density. 

• The 16 to 24 /im ratio shows significant variation with redshift. Most sources lie 
within the locus expected for starbursts, IR luminous galaxies, and AGN. A ratio 
> 1.4 appears likely to select predominantly sources at 1.1 < z < 1.6, where the 
the minimum between PAH emission peaks (as well as the 9.7 /im silicate absorption) 
redshifts into the 24 /im passband. About 5% of 16 /im sources meet this color selection, 
though few are selected as AGN by X-ray emission or power-law SEDs and so any AGN 
contribution must be heavily obscured. 

• Galaxy number counts show good agreement with previous surveys at similar wave- 
lengths (from ISOCAM and AKARI). The large number of sources and the two fields 
provide improvements in both the Poissonian errors and the effects of cosmic variance. 
AGN make substantial contribution to the number counts at bright fluxes. 

• We measure the integrated galaxy light at 16 /im, as a lower limit on the contribution 
to the monochromatic IR background, to be 2.2 ± 0.2 nW m~ 2 sr -1 . Extrapolating 
to fluxes 10 times fainter than the survey limit raises this value to 2.9 nW m -2 sr -1 , 
or 75% of the extragalactic background light. The contribution of sources which host 
AGN to the 16 /im EBL is ~ 15%, similar to the fraction found in 24 /im surveys. 

Future 15 /im results are expected from AKARI. Significantly more area has been cov- 
ered, which will allow improved study of rare sources. The Spitzer archive also contains 
many shallow (though small) fields which were obtained in parallel to spectroscopic obser- 
vations (Fajardo-Acosta et al. 2010, in prep.). Confusion-limited 16 /im observations were 
obtained by Spitzer in a small number of fields (e.g. Program ID=499; PI=Colbert). Finally, 
the recently launched WISE mission will include 12 /im images over the whole sky 

This work is based in part on observations made with the Spitzer Space Telescope, 
which is operated by the Jet Propulsion Laboratory, California Institute of Technology under 
NASA contract 1407. Support for this work was provided by NASA through an award 
issued by JPL/Caltech. This research has made use of the NASA/ IPAC Infrared Science 
Archive (IRS A), which is operated by the Jet Propulsion Laboratory, California Institute of 
Technology, under contract with the National Aeronautics and Space Administration. 
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Table I. GOODS-N, Spitzer data 



ID 


RA 


DEC 


Cov a 


fchl 


(? c hl 


fch2 


(J c h2 


fch3 


(T c / l 3 


fchA 


<J ch4 


/l6 


C"16 


/24 


<T24 




(deg) 


(deg) 




(MJy) 




(MJy) 


(>Jy) 


(MJy) 


(A*Jy) 


(MJy) 


(/iJy) 


OJy) 


(/*Jy) 


(/iJy) 


(^Jy) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(ii) 


(12) 


(13) 


(14) 


(15) 


(16) 





189.319641 


62.390766 


11 


84.80 


4.24 


60.40 


3.02 


55.20 


2.94 


40.20 


2.16 


145.1 


11.3 


225.0 


5.5 


1 


189.310349 


62.388977 


12 


11.50 


0.58 


14.50 


0.73 


19.00 


1.50 


12.90 


0.99 


55.3 


10.4 


200.0 


5.7 


2 


189.326126 


62.383057 


20 


38.30 


1.92 


49.10 


2.46 


55.30 


2.83 


39.30 


2.06 


107.2 


16.5 


483.0 


5.7 


3 


189.294266 


62.376274 


20 


40.00 


2.00 


47.50 


2.38 


39.60 


2.11 


44.20 


2.30 


240.4 


11.5 


386.0 


4.9 


4 


189.330811 


62.375191 


24 


10.10 


0.51 


6.98 


0.36 


7.48 


0.63 


6.27 


0.63 


53.9 


11.5 


26.6 


5.6 


5 


189.284683 


62.377537 


15 


33.50 


1.68 


30.20 


1.51 


21.00 


1.40 


20.50 


1.24 


96.0 


18.3 


89.1 


6.2 


6 


189.301178 


62.375061 


21 


10.90 


0.55 


10.50 


0.53 


8.67 


0.75 


11.60 


0.82 


75.5 


11.7 


60.5 


3.6 


7 


189.317520 


62.370590 


27 


28.40 


1.42 


24.10 


1.21 


20.70 


1.15 


34.00 


1.78 


129.3 


14.6 


133.0 


4.9 


8 


189.314499 


62.371925 


26 


1.16 


0.08 


1.13 


0.10 


0.47 


0.52 


3.13 


0.56 


56.7 


11.3 


50.7 


4.0 


9 


189.306412 


62.371181 


22 


50.70 


2.54 


37.50 


1.88 


32.20 


1.70 


30.10 


1.60 


339.4 


13.7 


436.0 


5.8 


10 


189.266174 


62.368641 


19 


7.67 


0.39 


5.57 


0.30 


4.81 


0.84 


4.93 


0.72 


62.8 


15.0 


18.3 


4.6 


11 


189.292343 


62.368526 


18 


32.10 


1.61 


29.90 


1.50 


21.80 


1.23 


23.50 


1.30 


130.8 


12.8 


125.0 


4.0 


12 


189.351440 


62.366024 


21 


242.00 


12.10 


199.00 


9.95 


131.00 


6.57 


214.00 


10.71 


365.8 


17.6 


395.0 


7.7 


13 


189.351089 


62.363888 


21 


18.80 


0.94 


14.00 


0.70 


12.40 


0.77 


14.00 


0.87 


62.8 


11.8 


60.5 


6.2 


14 


189.281219 


62.363293 


15 


47.40 


2.37 


51.00 


2.55 


52.10 


2.66 


59.60 


3.03 


129.2 


16.5 


240.0 


3.4 


15 


189.276138 


62.360157 


22 


101.00 


5.05 


103.00 


5.15 


117.00 


5.88 


123.00 


6.17 


219.1 


14.1 


303.0 


5.2 



Coverage in number of exposures (30 seconds integration per exposure). 



Table 2. GOODS-N, ancillary data 



Rcf. 1 X 2 SB flux 3 HBflux 3 S 4 Ratio 5 N c , 6 N c!16 7 Q 8 E 9 B 10 V 10 / 10 z 10 <r B w a v 10 a, 10 a z 10 Ids 11 

(10 -15 mW m -2 ) (mag) 



(17) 


(18) 


(19) 


(20) 


(21) 


(22) 


(23) 


(24) 


(25) 


(26) 


(27) 


(28) 


(29) 


(30) 


(31) 


(32) 


(33) 


(34) 


(35) 


(36) 


0.5820 


72 





0. 


.000 


0.000 





2.32 


1 


1 


1 





-999.00 


-999.00 


-999.00 


-999.00 


-999.00 


-999.00 


-999.00 


-999.00 





-99.0000 


-1 





0. 


.000 


0.000 





1.24 


1 


1 


3 





25.58 


26.24 


25.45 


24.75 


0.08 


0.28 


0.16 


0.09 


1 


-99.0000 


-1 








.000 


0.000 





3.09 


1 


1 


1 





26.85 


26.14 


25.80 


25.01 


1.45 


0.19 


0.17 


0.09 


1 


1.5220 


7.3 





0. 


.000 


0.000 





2.65 


2 


2 


5 





25.11 


24.43 


23.62 


23.09 


0.09 


0.08 


0.05 


0.03 


1 


0.6760 


72 





0. 


.000 


0.000 





1.76 


1 


1 


9 





24.71 


23.68 


22.82 


22.59 


0.04 


0.03 


0.02 


0.01 


1 


-99.0000 


-1 








.000 


0.000 





2.65 


1 


1 


1 





26.26 


24.98 


24.00 


23.00 


0.23 


0.11 


0.06 


0.03 


1 


-99.0000 


-1 





0. 


.000 


0.000 





2.04 


1 


1 


1 





27.90 


27.02 


25.50 


24.44 


0.50 


0.34 


0.11 


0.05 


1 


0.5045 


21 





0. 


.000 


0.000 





3.09 


1 


1 


1 





23.26 


22.19 


21.49 


21.17 


0.01 


0.01 


0.01 


0.01 


1 


-99.0000 


-1 








.000 


0.000 





3.72 


1 


1 


1 





25.82 


25.51 


25.04 


24.54 


0.08 


0.06 


0.05 


0.04 


1 


0.9746 


21 





0. 


.000 


0.000 





2.16 


1 


1 


1 





23.98 


23.14 


22.16 


21.65 


0.03 


0.02 


0.01 


0.01 


1 


1.0084 


21 








.000 


0.000 





3.61 


1 


1 


9 





24.38 


23.96 


23.27 


22.90 


0.02 


0.02 


0.01 


0.01 


1 


1.3644 


21 








.000 


0.000 





2.78 


1 


1 


1 





24.24 


23.82 


23.09 


22.45 


0.04 


0.03 


0.02 


0.01 


1 


0.2810 


21 





0. 


.000 


0.000 





3.27 


1 


2 


5 





20.70 


19.13 


18.39 


18.07 


0.01 


0.01 


0.01 


01.01 


1 


0.8615 


21 





0. 


.000 


0.000 





5.14 


2 


2 


7 





-999.00 


-999.00 


-999.00 


-999.00 


-999.00 


-999.00 


-999.00 


-9*g.OO 





1.4498 


21 


1 


8. 


.270 


14.120 





2.14 


1 


1 


1 





22.97 


22.77 


22.37 


21.98 


0.01 


0.01 


0.01 


o r oi 


1 


0.9032 


21 


1 


5. 


.760 


7.250 





3.32 


2 


1 


5 





22.61 


21.84 


21.31 


20.91 


0.01 


0.01 


0.01 


0.01 


1 



W-l no redshift; 10 -ICohen et al.1 i200($ ); ICohenl <200ll); 12 -iDawson etafl fcOOll); 21 -IWirth et al.l j2004h : 42 - IChapman et alj <2005h ISwinbank et al.1 [120041); 45 -IPope et al.1 
(200SftlMurphv et alj d200a> ,Spitzer IRS redshifts; 51 - iTreu et al. 1 1 |2005|) ; 61 - iReddv et al] [|2006h ; 72 — Stern et al. (in prep. ); 73 - iDaddi et alj feOOSl) ; 81 - iBarger ct al. I 1120081 ) 

2 Xray detection (l=detected, 0=not detected) 

3 Alexander et al. (2003) 

4 Star flag (l=star, 0=galaxy) 

5 Concentration (ratio of apflux in 6:2 pixel radii) 
6 Number of ch4 sources within 4 arcsec radius 
7 Number of catalog sources within 5.4 arcsec radius 

8 Quality Flag is assigned Bitwise: Bit = included in catalog; Bit 1 = more than one 16 micron source within 5.4 arcsec; possible confusion; Bit 2 = source concentration differs from 
that expected for point source; Bit 3 = S/N ratio < 5, but coverage indicates it should be higher 

9 Extcnd Flag (l=extended, 0=point source) 



Giavalisco et al. (2004) 

Number of i-band sources within 1 arcsec radius 
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Table 3. GOODS-S, Spitzer data 



ID 


RA 


DEC 


Cov a 


fchl 




fch2 




fch3 


^chZ 


fchA 


°"c/i4 


/l6 




/24 


CT24 




(deg) 


(deg) 




OJy) 


(^Jy) 


(A«Jy) 


(>Jy) 


OJy) 


(A*Jy) 


(MJy) 


(^Jy) 


(MJy) 


(MJy) 


(MJy) 


(^Jy) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(ii) 


(12) 


(13) 


(14) 


(15) 


(16) 





53.103943 


-27.663679 


2 


29.60 


1.48 


20.80 


1.05 


23.70 


1.37 


17.00 


1.10 


120.5 


23.7 


120.0 


5.7 


1 


53.142094 


-27.664965 


2 


33.70 


1.69 


26.00 


1.30 


21.40 


1.19 


22.30 


1.24 


252.7 


25.3 


161.0 


8.8 


2 


53.144711 


-27.666069 


3 


171.00 


8.55 


149.00 


7.45 


112.00 


5.62 


576.00 


28.81 


830.6 


39.8 


938.0 


11.5 


3 


53.110268 


-27.667072 


2 


26.90 


1.35 


22.30 


1.12 


16.30 


1.00 


16.80 


1.03 


208.5 


19.4 


143.0 


5.6 


4 


53.103607 


-27.666225 


4 


21.10 


1.06 


19.70 


0.99 


13.80 


0.92 


13.00 


0.90 


81.8 


14.6 


44.2 


4.9 


5 


53.102222 


-27.669676 


3 


46.90 


2.35 


35.90 


1.80 


33.60 


1.77 


40.40 


2.10 


251.0 


38.7 


322.0 


6.6 


6 


53.110622 


-27.669182 


3 


31.10 


1.56 


26.90 


1.35 


20.50 


1.17 


17.30 


1.04 


99.2 


17.9 


74.0 


5.4 


7 


53.123016 


-27.669071 


3 


102.00 


5.10 


63.70 


3.19 


60.70 


3.08 


40.60 


2.10 


201.5 


20.2 


256.0 


5.3 


8 


53.126549 


-27.669138 


3 


25.90 


1.30 


26.20 


1.31 


19.70 


1.11 


19.90 


1.14 


174.7 


15.9 


118.0 


5.1 


9 


53.129875 


-27.671673 


4 


61.40 


3.07 


43.80 


2.19 


45.70 


2.34 


39.40 


2.04 


199.7 


24.8 


203.0 


4.9 


10 


53.147671 


-27.672302 


4 


44.70 


2.24 


39.10 


1.96 


27.80 


1.47 


25.60 


1.38 


158.2 


24.8 


108.0 


4.9 


11 


53.102448 


-27.672508 


4 


28.80 


1.44 


22.20 


1.11 


16.80 


0.99 


16.30 


0.98 


169.0 


16.8 


162.0 


5.5 


12 


53.109512 


-27.674276 


4 


28.50 


1.43 


22.10 


1.11 


19.20 


1.08 


22.10 


1.23 


180.3 


17.4 


175.0 


11.9 


13 


53.114948 


-27.674074 


4 


25.50 


1.28 


18.20 


0.91 


13.90 


0.85 


12.10 


0.80 


107.8 


17.0 


67.8 


5.4 


14 


53.127071 


-27.675257 


4 


33.80 


1.69 


32.10 


1.61 


25.40 


1.35 


19.10 


1.09 


116.0 


16.5 


57.1 


3.9 


15 


53.090645 


-27.675821 


3 


42.60 


2.13 


41.00 


2.05 


33.30 


1.74 


84.40 


4.25 


215.6 


19.8 


272.0 


4.0 



( a ) Coverage in number of exposures (60 seconds per exposure). 



Table 4. GOODS-S, ancillary data 



Ids 1 



2spcc Ref. 1 X 2 SB flux 3 HBflux 3 S 4 Ratio 5 N ci 6 N cil6 7 Q 8 E 9 B 10 V 10 I 10 z 10 a B 10 a v 10 cr 7 10 a z 10 

(10- 15 mW m~ 2 ) (mag) 
(17) (18) (19) (20) (21) (22) (23) (24) (25) (26) (27) (28) (29) (30) (31) (32) (33) (34) (35) (36) 



0.6658 


12 





0. 


.000 


0.000 





2.53 


1 


1 


17 





25.32 


23.39 


22.42 


22.09 


0.01 


0.02 


0.02 


0. 


.01 


1 


1.1120 


30 








.000 


0.000 





3.31 


2 


1 


27 





23.07 


22.54 


21.96 


21.52 


0.01 


0.01 


0.01 





.01 


1 


0.2140 


42 








.000 


0.000 





2.72 


1 


1 


1 





20.50 


19.34 


18.79 


18.57 


0.01 


0.01 


0.01 


0. 


.01 


2 


0.6243 


1 





0. 


.000 


0.000 





2.49 


1 


1 


17 





26.43 


24.76 


23.77 


23.06 


0.01 


0.05 


0.04 


0. 


.02 


1 


-99.0000 


-1 








.000 


0.000 





2.70 


1 


1 


1 





26.82 


23.40 


25.00 


24.27 


0.01 


0.02 


0.13 





.08 


1 


-99.0000 


-1 


1 





.000 


1.300 





2.58 


2 


1 


3 





24.53 


23.66 


22.55 


21.91 


0.04 


0.02 


0.01 





.01 


1 


-99.0000 


-1 





0. 


.000 


0.000 





3.21 


1 


2 


3 





-999.00 


27.25 


27.28 


25.01 


-999.00 


0.43 


0.76 


0. 


.10 


2 


0.7343 


12 





0. 


.000 


0.000 





2.30 


1 


1 


1 





23.46 


22.25 


21.08 


20.64 


0.03 


0.01 


0.01 


0. 


.01 


1 


1.3570 


30 








.000 


0.000 





1.96 


1 


1 


1 





24.55 


24.08 


23.30 


22.60 


0.05 


0.03 


0.03 


0. 


.02 


1 


0.5659 


12 





0. 


.000 


0.000 





2.53 


1 


1 


1 





22.93 


21.69 


20.86 


20.52 


0.02 


0.01 


0.01 


0. 


.01 


1 


1.1240 


30 








.000 


0.000 





2.15 


1 


1 


1 





26.64 


25.34 


24.11 


23.01 


0.28 


0.07 


0.04 


0. 


.02 


1 


1.0368 


42 








.000 


0.000 





2.88 


1 


1 


1 





25.21 


24.09 


23.14 


22.48 


0.09 


0.02 


0.02 


0. 


.01 


1 


1.0387 


42 





0. 


.000 


0.000 





4.53 


1 


1 


5 





23.33 


22.93 


22.34 


22.03 


0.02 


0.01 


0.01 


0. 


.01 


1 


-99.0000 


-1 





0. 


.000 


0.000 





2.15 


1 


1 


1 





23.68 


23.23 


22.46 


21.92 


0.04 


0.02 


0.02 


0. 


.01 


1 


-99.0000 


-1 








.000 


0.000 





2.43 


1 


1 


1 





24.69 


24.16 


23.37 


22.64 


0.08 


0.05 


0.04 





.02 


1 


0.4260 


72 








.000 


0.000 





2.85 


1 


1 


1 





22.06 


20.98 


20.47 


20.20 


0.01 


0.01 


0.01 





.01 


1 



CO 
CO 



x l = LeFevre et al. (2004) ; 2 - Szokoly et al. (2004) ; 3 - Croom et al. (2001) ; 9 - K20 Survey, Mignoli et al. (2005) ; 11 - Strolger at al. (2004) ; 30 - Vanzella et 
al. (20 08) ; 41,42 - IPopesso et ail fcOOgl l lBalestra et ail ll2010h ; 62 - Doherty et al. (2005) ; 70 - Ravikumar et al. (2007) ; 72 - Stern et al. (in prep.) ; 81 - iKriek et ail 
{20081 ) 

2 Xray detection (l=detected, 0=not detected) 

3 Luo ct al. (2003) 

4 Star flag (l=star, 0=galaxy) 

5 Concentration (ratio of apflux in 6:2 pixel radii) 
6 Number of ch4 sources within 4 arcsec radius 
7 Number of catalog sources within 5.4 arcsec radius 

8 Quality Flag is assigned Bitwise: Bit = included in catalog; Bit 1 = more than one 16 micron source within 5.4 arcsec; possible confusion; Bit 2 = source concentration 
differs from that expected for point source; Bit 3 = S/N ratio < 5, but coverage indicates it should be higher; Bit 4 = source with coverage of two exposures instead of 
three 



9 Extend Flag (l=extended, O=point source) 
10 Giavalisco ct al. (2004) 

11 Numbcr of i-band sources within 1 arcscc radius 
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Table 5. Numbers of Sources 



Field 


Sources 


Hard Xray 


other power-law AGN 


Blue (/ 16 //24 > 1.4) a 


North 


840 


96 


9 


55 (25) 


South 


476 


58 


14 


51 (34) 



a The number of sources with /16//24 > 1-4 and a 5cr detection at 24 /zm, with either no 
spectroscopic redshift or a redshift 1.1 < z < 1.6, indicating possible silicate absorption. 
The number in parenthesis indicates how many of these sources have a quality flag value of 
1. 
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Table 6. Sensitivities 



Quantity 


North 


South 


UDF 


Min Itime per pixel (s) a 


300 


120 


600 


Planned Itime per pixel (s) 


600 


120 


960 


Mode Itime per pixel (s) 


630 


270 


2100 


Min predicted depth b (/xJy) 


80 


125 


55 


Mode predicted depth (/xJy) 


48 


80 


30 


Mode measured Depth c (/iJy) 


45 


80 


30 



a Minimum exposure time per pixel for objects in the catalog 

b Minimum predicted 5a depth using the exposure time calculator on the SSC website 
c Mode of measured 5a depth 
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Table 7. GOODS-North Galaxy Number Counts 



Slow 


Shigh 


Savg 


■N galaxies 


Observed dN/dS 


Corrected dN/dS 


Corrected 5 dN/dS 


(mJy) 


(mJy) 


(mJy) 




(deg- 2 mJy- 1 ) 


(deg- 2 mJy- 1 ) 


(deg~ 2 mJy -1 ) 


0.025 


0.040 


0.034 


13 








0.040 


0.056 


0.048 


106 


3.0xl0 5 


4.5xl0 5 


1.6xl0 4 


0.080 


0.112 


0.096 


151 


1.1x10 s 


1.3xl0 5 


9.6xl0 3 


0.112 


0.159 


0.136 


128 


6.6xl0 4 


7.1xl0 4 


6.0xl0 3 


0.159 


0.224 


0.192 


93 


3.4xl0 4 


3.4xl0 4 


3.5xl0 3 


0.224 


0.317 


0.271 


61 


1.6xl0 4 


1.6xl0 4 


2.0xl0 3 


0.317 


0.448 


0.382 


45 


8.2xl0 3 


8.2xl0 3 


1.4xl0 3 


0.448 


0.632 


0.540 


22 


2.8 xlO 3 


2.8xl0 3 


7.4xl0 2 


0.632 


0.893 


0.763 


11 


l.OxlO 3 


l.OxlO 3 


4.0xl0 2 


0.893 


1.262 


1.078 


3 


1.9xl0 2 


1.9xl0 2 


1.6xl0 2 


1.262 


1.783 


1.522 


2 


9.1X10 1 


9.1X10 1 


7.3X10 1 


1.783 


2.518 


2.150 











2.518 


3.557 


3.037 
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Table 8. GOODS-South Galaxy Number Counts 



Slow 

(mJy) 


Shigh 

(mJy) 


<7 

u avg 

(mJy) 


■N galaxies 


Observed dN/dS 
(deg~ 2 mJy -1 ) 


Corrected dN/dS 
(deg~ 2 mJy -1 ) 


Corrected^ dN/dS 
(deg -2 mJy -1 ) 


0.025 


0.040 


0.034 


14 


3.1xl0 5 


6.2xl0 5 


1.4xl0 5 


0.040 


0.056 


0.048 


20 


2.6xl0 5 


3.7xl0 5 


8.4xl0 4 


0.056 


0.080 


0.068 


38 


1.5xl0 5 


1.8xl0 5 


2.8xl0 4 


0.080 


0.112 


0.096 


107 


9.2xl0 4 


1.2xl0 5 


l.OxlO 4 


0.112 


0.159 


0.136 


109 


5.9xl0 4 


6.6xl0 4 


6.0 xlO 3 


0.159 


0.224 


0.192 


83 


3.2xl0 4 


3.4xl0 4 


3.6xl0 3 


0.224 


0.317 


0.271 


45 


1.2xl0 4 


1.2xl0 4 


2.1xl0 3 


0.317 


0.448 


0.382 


21 


4.0xl0 3 


4.0xl0 3 


l.lxlO 3 


0.448 


0.632 


0.540 


11 


1.5xl0 3 


1.5xl0 3 


6.0xl0 2 


0.632 


0.893 


0.763 


10 


9.6xl0 2 


9.6xl0 2 


4.1xl0 2 


0.893 


1.262 


1.078 


3 


2.0xl0 2 


2.0xl0 2 


1.6xl0 2 


1.262 


1.783 


1.522 


1 


4.8X10 1 


4.8X10 1 


3.9X10 1 


1.783 


2.518 


2.150 


2 


6.8X10 1 


6.8X10 1 


5.5X10 1 


2.518 


3.557 


3.037 


1 


2.4X10 1 


2.4X10 1 


1.9X10 1 
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a (2000) 



1,1 189.0 188.9 



53.2 53.1 
a (2000) 



Fig. 1. — We show the mosaic of the North (left) and South (right) GOODS fields. A 
smaller region is shown below to give a better indication of crowding and image quality; for 
the South, the region includes the UDF, which has greater sensitivity than the rest of the 
field. 



4 5 6 7 8 9 10 20 
Wavelength (Jim) 



Fig. 2. — Spitzer photometry for 15 stars in GOODS (circles), normalized to the the average 
of the 6-band photometry. Data are compared to IRS spectra of stars with stellar types A 
through M (grey spectra) and the iKuruca ( 19791 ) model for an AOV star (dashed spectrum). 
The IRS spectra are taken from the SASS survey (Ardila et al. in prep). Error bars are 
small for the IRAC data points; typical error bars for the 16 and 24 /im photometry are 
shown next to the data. 



-41 - 




Fig. 3. — The flux density measured with Spitzer PUI in th e 16 fim bandpass compared to 
that measured with ISOCAM in the 15 /zm bandpass (LW3 lAussel et al.lll999l . their Table 
3) for objects detected by both instruments. The dot-dashed line indicates equal flux, and 
the solid line indicates the best fit offset of a factor of 1.13. 
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Fig. 4. — The ratio of the Spitzer 16 fim to ISOCAM 15 /im (LW3) flux densities as a 
function of spectroscopic redshift. The four brightest sources are circled. Also plotted are 
predicted ratios based on Spitzer IRS spectra of templa te galaxies: the extreme silicate- 
absorption galaxy IRAS F00183-7111 dSpoon et alibooi solid line), UGC5101, a ULIRG 



with con siderable 9. 7/xm absor ption (lArmus et al.ll2004l . dashed line), the prototy pical AGN 



Mrk2 31 (jWeedman et al.ll2005l dotted line), the typical quasar PG1501 + 106 from (jHao et al. 



20051 . triple-dot-dash ed line) and the ave rage mid-IR SED of all starburst galaxies in the IRS 



GTO program from ( Brandl et al. 20061 . dot-dashed line). 
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Fig. 5. — The flux density measured with Spitzer PUI in the 16 lira bandpass c ompared to 
that measured with AKARI IRC in the 15 /mi bandpass ( jBurgarella et al.l 120091 ) for objects 
detected by both instruments. The dot-dashed line indicates equal flux, and the solid line 
indicates the best-fit offset which corresponds to a factor of 1.3. 
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200 




Predicted 5 sigma limit (|oJy) 



Fig. 6. — The area covered at expected sensitivity limits (in /iJy) in the North (solid line) 
and South (dashed line). Sensitivity predictions are based upon the SSC exposure time 
calculator (see text). 
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Redshift Redshif'l 



Fig. 7. — The distribution of spectroscopic redshifts (filled histogram) for sources associated 
with 16 /iin detections in the North (left) and South (right), compared to the distribution 
of sources associated with 24 /xm detections (dashed line) and the full spectroscopic sample 
(solid line; see text) 
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Fig. 8. — The flux density of 16 /zm detected sources versus spectroscopic redshift for GOODS 
North (black circles) and South (grey squares). 



-47- 



\ i i i i r 



0.4 



# AGN (xray) 

• AGN (xray or powerlaw) 



o 

•a 

o 



5 

< 



0.3 



0.2 



0.1 



0.0 



it 



it 



l * 



J A I I I I L 



10 



100 

16 |im flux density (|iJy) 



1000 



Fig. 9. — The fraction of 16 /xm sources identified as AGN (filled circles) by the detection of 
hard X-rays (open star symbols) and/or power law slope. 
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Fig. 10. — The ratio of 16 to 24 /im flux densities as a function of spectroscopic redshift 
for high quality (flag=l; black points) and all sources (grey points). The shaded region 
indi cates the range of values expected from the starburst templates of Brandl et al. (2006) 



and ISmith et al.l (120071 ) . The ratio for Arp220 (using the spectrum from lArmus et al.ll2007l ) 



is plotted (solid line), as is the value for a power law spectrum with f v ~ v~ 2 (dashed line). 
Uncertainties are plotted for objects with colors more than la bluer than expected from the 
local starburst templates. 
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Fig. 11. — As in Figure [101 we show the flux ratio for high quality sources (flag=l; black 
points) and highlight objects identified as possible AGN through hard X-ray detection or 
power law slope (star-shaped symbols, see section |4~T]) . 
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Fig. 12. — Histogram of 16 to 24 /zm flux ratios for all objects (solid line) and objects without 
spectroscopic redshift (left; filled histogram for objects with photometry quality flag equal 
to 1; dashed line for all objects without redshift) and AGN (right; dashed line). 
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Fig. 13. — Differential 16 fim number counts measured by Spitzer (filled red symbols have 
completeness correction applied, open red symbols do not), AKARI (blue symbols), and ISO 
(dark grey squares). The Euclidean slope has been removed. Poissonian error bars are 



shown. ISOCAM points include HDF-N, HDF-S, and the Marano su rveys (lElbaz et al. 



and the references therein); the gravitational lensing cluster survey ( lAltieri et al.lll999l ); the 



1999 



European Large- Area ISO Surve y (ELIAS-Sl;lGruppioni et al.ll2002l ); and the Lockman Deep 
and Lockman Shallow surveys ( Rodighiero et alT2004 ). We plot AKARI poi nts from the 



completeness-corrected co unts measured in the N o rth ec liptic pole (filled squares; iWada et al. 



20071 : IPearson et al.ll2010f ). The Burgarella et all (12009( 1 survey of GOODS-South (filled blue 
triangles are counts with completeness correction, open blue triangles are the counts without) 
and we exclude bins with less than 20% completeness. We also show the com pleteness- 
corrected AKARI counts with lensed objects in Abell 2218 (IHopwood et al.ll2010l ). We also 
show the fit to the combined Spitzer 16 fim counts from the two GOODS fields (dotted line; 
see Section |4~4I) 
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Fig. 14. — Differential 16 /zm number counts for GOODS North and South combined (filled 
circles). We show the contribution of sources identified as AGN (star-shaped symbols). We 
also show the contribution of objects with 16 to 24 /iin flux ratios > 1.4, and spectroscopic 
redshifts 1 < z < 2, indicating possible silicate absorption (open circles). 
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Fig. 15. — The integrated galaxy light (IGL) as a function of sensitivity limit. The IGL 
is extrapolated below the UDF sensitivity (dotted vertical line) assuming constant faint-end 
slope. The uncertainty (dotted curve) is estimated by calculating the IGL fo r the la upper 



and l ower limits to the number counts. The upper limit to the IGL is shown (jRenault et al. 



20011 ) as the solid horizontal line. 



